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Abstract. If the Galactic WMAP radio haze, as recently confirmed by Planck, is produced 
by dark matter annihilation or decay, similar diffuse radio halos should exist around other 
galaxies with physical properties comparable to the Milky Way. If instead the haze is due 
to an astrophysical mechanism peculiar to the Milky Way or to a transient event, a similar 
halo need not exist around all Milky Way "twins". We use radio observations of 66 spiral 
galaxies to test the dark matter origin of the haze. We select galaxies based on morphological 
type and maximal rotational velocity, and obtain their luminosities from a 1.49 GHz catalog 
and additional radio observations at other frequencies. We find many instances of galaxies 
with radio emission that is less than 5% as bright as naively expected from dark matter 
models that could produce the Milky Way haze, and at least 3 galaxies that are less than 
1% as bright as expected, assuming dark matter distributions, magnetic fields, and cosmic 
ray propagation parameters equal to those of the Milky Way. For reasonable ranges for the 
variation of these parameters, we estimate the fraction of galaxies that should be expected 
to be significantly less bright in radio, and argue that this is marginally compatible with the 
observed distribution. While our findings therefore cannot rule out a dark matter origin for 
the radio haze at this time, we find numerous examples (including the Andromeda Galaxy) 
where, if dark matter is indeed the origin of the Milky Way haze, some mechanism must be 
in place to suppress the corresponding haze of the external galaxy. We point out that Planck 
data will offer opportunities to improve this type of constraint in a highly relevant frequency 
range and for a potentially larger set of candidate galaxies. 
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1 Introduction 

Data from NASA's Wilkinson Microwave Anistropy Probe (WMAP) have fostered great ad- 
vances in both precision cosmology and in the understanding of processes associated with 
the interstellar medium of our own Galaxy. WMAP also brought about unexpected observa- 
tional results, including the detection of spinning dust emission [1-4] and a mysterious diffuse 
microwave emission extending 20° around the Galactic Center, known as the "WMAP haze" 
[5]. While originally associated with free- free emission [5], follow-up analyses showed that 
the haze spectrum was too soft to be attributed to free-free emission, but too hard to match 
the diffuse Galactic synchrotron emission produced by relativistic cosmic ray electrons [6]. 

While some authors have criticized the original analysis procedure [7] , or the synchrotron 
nature of the haze [8], both arguments have been rebutted in detail [9, 10]. Excitingly, the 
Planck mission has recently presented conclusive evidence in favor of the existence of the mi- 
crowave haze [11], employing superior background rejection methods to those available to the 
WMAP team. Specifically, the Planck data analysis proceeded along two distinct component 
separation techniques, one following the original WMAP template analysis of Ref. [6], and a 
more sophisticated Bayesian approach based on Gibbs sampling. Both techniques confirmed 
the existence of an anomalous diffuse haze in the Galactic Center (GC) region across mul- 
tiple frequencies, with a spectral index (3h = —2.55 it 0.05, consistent with that found from 
WMAP data. Furthermore, future Planck studies will analyze the polarization of the haze, 
providing an even more stringent test of the synchrotron nature of this emission. 

Several authors have considered dark matter annihilation as a possible source of the 
haze (see e.g. Ref. [12-17]). The dark matter interpretation is consistent with theoretically 
well motivated ranges for the mass and annihilation cross section of weakly interacting mas- 
sive particles, common to many extensions of the standard model of particle physics [18-20]. 
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Furthermore, it was noted recently that the WMAP haze is consistent with the annihila- 
tions of light (m ~ 10 GeV) dark matter candidates which are also capable of accounting 
for a number of other reported observations [16], including direct detection signals from 
DAMA/LIBRA [21], CoGeNT [22, 23] and CRESST-II [24], as well as the spatially extended 
spherically symmetric excess of 7-rays observed around the GC by the Fermi-LAT [25-27] , 
the extremely hard synchrotron spectrum observed in multiple filamentary arcs surrounding 
the GC [28] , and the isotropic radio excess observed by the ARCADE collaboration [29] (see 
Ref. [30] for further discussion of these observations and experiments). 

Recent observations indicate that the observed morphology of the microwave excess 
haze exhibits a strong correlation with a diffuse feature discovered at GeV energies by the 
Fermi Large Area Telescope, known as the Fermi bubbles [31] (previously known as "Fermi 
haze" [32], see also [33]). This suggests a common origin for the observed emissions at radio 
and gamma-ray energies [34]. Hypotheses as to the origin of this additional cosmic ray 
population include anomalous enhancements to the supernova activity in the relevant region 
and time [35], peculiar cosmic ray propagation in the region of interest, perhaps associated 
with Galactic winds [36], central nuclear activity fueling cosmic ray jets [37], or an additional 
population of cosmic ray electrons distinct from the ordinary diffuse cosmic ray population 
in the Galactic plane [14]. High mass dark matter models may also simultaneously explain 
both the features of the WMAP haze and Fermi bubbles, although highly anisotropic diffusion 
mechanisms are needed to produce the observed emission morphology [15]. 

In this paper, we point out that if the WMAP haze results from dark matter annihila- 
tion, a radio signal with a similar luminosity and morphology should be emitted from other 
galaxies, so long as these galaxies are comparable to our own in terms of size and other 
physical characteristics. This comparison may not hold in scenarios where the haze is of 
astrophysical origin (e.g. in the case that the haze is powered by episodic nuclear activity 
in the GC). The purpose of this study is to identify a suitable set of galaxies similar to 
the Milky Way, estimate the radio "haze" expected from these candidate galaxies assuming 
a dark matter annihilation origin and benchmark values for the parameters governing the 
propagation of charged leptons in each galaxy, and then employ radio observations to obtain 
constraints on the dark matter annihilation origin of the Milky Way haze. 

In order to interpret the results of this comparison, it is necessary to estimate the 
uncertainties in the theoretical expectation for the dark matter radio haze in external galaxies. 
Relatively little is known about propagation of cosmic rays, magnetic fields, or dark matter 
distribution in external galaxies. Any or all of these characteristics could vary considerably 
from galaxy-to-galaxy, even among those galaxies that are similar in size and morphology 
to the Milky Way. We therefore study in detail the expected range of variation in the radio 
luminosity of the dark matter radio haze originating from such variations. 

Our study is structured as follows. In the next section, we describe the selection criteria 
we employ to select the relevant set of galaxies to investigate. In Sec. 3, we address the 
systematic uncertainties on the estimate of the haze luminosity. In Sec. 4 we produce a Monte 
Carlo model indicating the expected variation in the dark matter induced synchrotron signal 
due to these systematic uncertainties. In Sec. 5, we compare our predictions with radio data. 
Lastly, we discuss implications of our results and future directions in Sec. 6. 
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Type 


t 


i 

[deg.] 


frot 

[km / s] 


S 
[mJy] 


Sp 

[mJy] 


Dist. 

[Mpc] 


Rank^ 


Rank2kpc 


Avg. Rank 


NGC 4448 


SBab 


1.8±0.6 


69.00 [39] 


221.54 


1 


0.9 


13.1 


1 


2 


1.5 


NGC 4698 


Sab 


1.7±1.0 


73.44 


201.13 


0.6 




21.9 


2 




2 


NGC 4394 


SBb 


3 ±0.4 


16.55 


255.13 


0.7 


0.6 


21.9 


3 


3 


3 


NGC 7814 


Sab 


2±0.2 


90.00 [40] 


230.9 


1.1 


1.1 


25.8 


4 


6 


5 


NGC 1350 


Sab 


1.9±0.6 


64.79 


199.87 


1.1 


0.8 


30.2 


5 


7 


6 


NGC 0224 


Sb 


3±0.4 


72.17 


256.7 


8400 


14.5 


0.7 


12 


1 


6.5 


NGC 2683 


Sb 


3±0.3 


82.79 


202.92 


65.9 


14.6 


8 


13 


5 


9 



Table 1. Candidate galaxies from the Condon Atlas, t, i , v^ot , S, Sp, refer to the de Vaucouleurs 
morphological type, inclination, maximum rotational velocity, total flux, and peak flux, respectively. 
Ranki (Rank2kpc) ranks the total (respectively, 2 kpc off disk) luminosity, in ascending order. Unless 
otherwise noted, inclination, type, and rotational velocity are from the HyperLeda database. 

2 Galaxy Sample Selection 

In order to identify galaxies with size and morphology similar to the Milky Way, we utilize 
the 1.49 GHz radio data from the Condon Atlas of Spiral Galaxies (hereafter refered to as 
CA) [38]. We obtain morphological information for each galaxy by cross-referencing the CA 
with the HyperLeda database^, and restrict our analysis to galaxies with a morphological 
type 1.5 < t < 4.0 ^, which limits the uncertainties associated with edge-on galaxies, for 
which morphological properties can be difficult to determine. We note that this conservative 
cut loosely includes type Sab-Sbc galaxies. 

Apart from well-studied galaxies such as Andromeda, no specific information on the dark 
matter density profile is available for the galaxies in the CA. Most of these galaxies, however, 
have estimates for the values of the maximum rotational velocity, given by observations of 
21 cm line widths. These values are highly correlated with the total enclosed mass in each 
galaxy. We place a cut on the maximum rotational velocity of 180 km/s< Vrot < 280 km/s 
to ensure a mass and size comparible with the Milky Way. 

After each of the above cuts is applied, the resulting sample contains 66 galaxies. We 
then rank the galaxy candidates based on the observed radio flux in two ways. First, we 
calculate the total radio luminosity at 1.49 GHz given the total observed radio flux (S) 
reported in the CA and the best-fit distance estimate as given by the NASA/IPAC Extra- 
galactic Database (NED)^. Second, we calculate the luminosity 2 kpc above or below the 
Galactic plane by using the peak flux measurement (Sp) reported by the CA, as well as the 
elliptical Gaussian defined by the FWHM at 1.49 Ghz along the major and minor axes, with 
the major axis assumed to lie along the Galactic plane as calculated by optical data. The 
resulting list is then sorted in ascending order of luminosity, and the two rankings are then 
averaged. In Table 1, we list the seven top-ranking galaxies, which we will use throughout 
this analysis. Notably, this list of selected galaxies includes Andromeda (NGC 0244), which 
has well studied physical properties which identify it as a veritable Milky Way "twin" . 

In order to estimate the physical size of the cosmic ray diffusion region, we use as a 
proxy the projected size of the major axis of each galaxy, based on the distances obtained 

^The HyperLeda database may be found at http://leda.univ-lyonl.fr/ 

■^Mapping between Hubble classification and de Vaucouleurs morphological type index may be found at 
http:/ /leda.univ-lyonl.fr/leda/param/t.html 

^The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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Name NED Distance Major Axis Axes Ratio 
[Mpc] [kpc] 



NGC 0224 
NGC 1350 
NGC 2683 
NGC 4394 
NGC 4448 
NGC 4698 
NGC 7814 



0.7 
20.9 
10.2 
16.8 
13.0 
23.7 
17.2 



60.2 
52.4 
46.9 
28.4 
10.0 
43.8 
36.4 



2.55 ± 0.29 
1.99 ± 0.15 
3.51 ± 0.42 
1.04 ± 0.07 
1.53 ± 0.15 
2.45 ± 0.16 
2.32 ± 0.14 



Table 2. Distances and sizes of the candidate sample galaxies. 



from NED and an angular size from HyperLeda, defined as the isophotal level 25 mag/arcsec^ 
in the B-band [41]. In Table 2, we list the mean distance for the galaxy sample, as given in 
the NED database. For each galaxy in the sample, we extrapolate the physical size of the 
Galaxy's major axes. We also report the ratio of the major-to-minor axis in the last column, 
although we assume cylindrical symmetry for the simulated diffusion region. 

We conclude that all of our candidate galaxies have comparable diffusion region sizes, 
with the possible exceptions of NGC 4448 and NGC 4394, which exhibit a potentially smaller 
radius and a smaller major-to-minor axis ratio. These galaxies also have the largest distance 
uncertainties, directly affecting the reported semi-major axis. 

3 The Dark Matter Haze in External Spiral Galaxies: Systematic Uncer- 



Although we employ morphological information to choose a sample of galaxies similar to the 
Milky Way, there are several other parameters capable of impacting the radio luminosity of 
a dark matter haze which can vary significantly from galaxy-to-galaxy. Generically, these 
parameters break down into three categories: 

1. parameters that control the dark matter density profile, 

2. parameters that control the diffusion of cosmic rays, and 

3. parameters that control the magnetic field strength in each sample galaxy. 

In Table 3 we show the default, maximum and minimum value for each parameter used 
throughout this work, and in Figure 2 we show the effect of changes in each individual 
parameter on the synchrotron luminosity due to dark matter annihilation at 1.49 GHz. In 
the following, we discuss the variations in the radio haze produced by of these parameters in 
detail. 

3.1 Dark Matter Density Profile 

As the dark matter annihilation rate is proportional to the square of the dark matter density, 
the total synchrotron luminosity depends not only on the total dark matter mas (which 
can be reasonably estimated from galactic rotation measurements) but also on the density 
distribution of dark matter in each galaxy. In this work, we assume a generalized Navarro- 
Frenk- White (NFW) profile [42] with a radial density profile given by: 
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2 
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5 


0.30 


0.06 
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2 


22 


11 


44 


a 




1.5 


1.0 


0.67 


1.5 



Table 3. Nuisance parameters employed to estimate the range for a dark matter induced haze in 
external spiral galaxies. Note that we also place an additional constraint on each simulation that the 
diffusion height must not exceed the simulated diffusion radius. 

where is a scale radius which governs the turnover to the r~" profile, ps is a dark matter 
density normalization constant, and a governs the inner slope of the dark matter density. In 
the original NFW profile, a = 1, in agreement with the results of certain dark matter-only 
simulations of structure formation, e.g. Ref. [43]. However, recent simulations which include 
the effects of baryons on galaxy evolution have found that the cooling and contraction of 
the baryonic density profile tend to steepen the dark matter density profile in the central 
regions [44]. This effect, known as baryonic contraction, generally leads to values of a in a 
range between 1.2-1.7, depending on the degree of contraction. In contrast, other numerical 
simulations which include a large degree of baryonic feedback have found that baryons can 
flatten the inner slopes of the dark matter profile, producing distributions with a flat density 
core as large as 1 kpc in size [45]. In our simulations we allow for r^ to change by a factor 
of 2 from our default choice of 22.0 kpc, and we allow a to change by a factor of 1.5 from a 
default value of 1.0. We normalize the dark matter density profile the dark matter density 
at the solar position (R0 = 8.5 kpc) to a range within a factor of 5 of the central value of 
0.3 GeV cm-3 [46]. 

Figure 1 (top) illustrates that the haze luminosity depends strongly on the assumed 
value for a, increasing precipitously for profiles with q > 1. We note that this will affect our 
sample producing a logarithmically asymmetric distribution in luminosity and favoring large 
radio fluxes. We also note only a marginal dependence on the value of as the annihilation 
which produces the radio haze is largely confined to within the inner several kpc around the 
center of each galaxy. 

3.2 Cosmic Ray Diffusion 

As charged particles propagate through the galaxy, they are deflected by the turbulent mag- 
netic field structure and proceed in a random walk that can be described as a standard 
diffusive process. Since Galactic magnetic fields are thought to be produced primarily by 
magnetohydrodynamic turbulence, they are highly dependent on astrophysical parameters 
such as the distribution of molecular clouds, ionized gas, star formation rates, etc. which 
vary wildly between galaxies [47]. Additionally, the nature of charged particle propagation 
depends sensitively on the scale of the magnetogydronamic fluctuations compared to the 
gyroradius of the charged particle in the magnetic field. This is a local phenomena which is 
not uniform across different regions of the Galaxy. 
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Figure 1. The effect of variations in the nuisance parameters on the luminosity of the dark matter 
synchrotron haze at 1.49 GHz in units of the benchmark value, Lmw- See Table 3 and text for more 
details. 

Although an exact solution to the diffusion equation for a galaxy is not attainable, cos- 
mic ray propagation codes exist that allow for realistic estimates of charged particle diffusion 
in the Milky Way. The predictions of these codes can be tested by observable quantities such 
as the ratios of radioactive secondaries to primary species observed at the solar position, 
large scale gamma-ray emission, and other diffuse Galactic radiation backgrounds. The Gal- 
prop code stands at the forefront of this effort, performing numerical calculations of particle 
propagation and of the resulting multi-wavelength emission, taking into account effects such 
as the diffusion, reacceleration, and convection of cosmic rays, the intensity as well as the 
morphology of the Galactic magnetic and interstellar radiation fields, and the morphology of 
molecular gas [48]. In this analysis, we employ Galprop (v. 54) to study the effect of varying 
four parameters, which previous analyses of the WMAP haze in the Milky Way determined 
to have the largest effect on the synchrotron signal due to dark matter annihilation [17]. 
Specifically we examine changes in the mean diffusion constant, Dq, whose benchmark value 
we take to be 1.0 x 10^^ cm^ s~^, the radius (-Rdiff) and half-height (/idiff) of the diffusion 
region, which we take to be 20 kpc and 16 kpc for our central values, and the normalization 
of the inter-stellar radiation field, which we take to be unity in units of the default Galprop 
model. 
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Figure 1 (middle) examines the changes in the mean Milky Way luminosity for variations 
in our diffusion setup, finding only an minimal effect from changes in the size of the diffusion 
region, on the order of 10% for regions double the default size. However, changes in the 
diffusion constant produce a larger effect, especially in the regime of low diffusion constants 
where electrons created by dark matter annihilations are efficiently trapped in the galaxy. 
Large increases in the strength of the interstellar radiation field also greatly decrease the 
luminosity of the synchrotron haze, as charged leptons lose a significant fraction of their 
power to inverse Compton scattering as opposed to synchrotron. 

3.3 Galactic Magnetic Fields 

In addition to influencing the propagation of charged cosmic rays. Galactic magnetic fields 
also control the intensity and spectrum of synchrotron radiation produced by cosmic ray 
electrons as they travel through the interstellar medium. Unfortunately, the mean intensity 
of magnetic fields is highly uncertain, even in the case of the Milky Way. However, constraints 
from the rotation measures of pulsars place the mean magnetic field strength in the Milky Way 
to be approximately 4-6 nG [49] . In this paper, we adopt the a magnetic field model described 
by the functional form B{r, z) = Bq exp(— r/rg) exp(— z/zq), and vary the parameters of this 
expression as described in Table 3. 

In Figure 1 (bottom) we show the variation in the dark matter synchrotron luminosity 
for changes in the strength of the magnetic field in a given galaxy. We find only moderate 
changes when we shift the overall amplitude of the magnetic field, which is primarily due 
to the large fraction of electron energy which is lost to synchrotron for all magnetic field 
strengths. The dependence of the total synchrotron intensity on the radial and zenith scale 
of the magnetic field is even more negligible, due to the fact that the synchrotron intensity 
is dominated by the very inner region around each galaxies' center. 

4 Monte Carlo Modeling of Systematic Uncertainties 

In order to estimate a probability distribution for the luminosity of each external galaxy, we 
produce a Monte Carlo sample, varying each parameter listed above within the minimum 
and maximum values provided in Table 3. In the case where the variation is less than or 
equal to a factor of two alteration in the parameter, we sample the parameter space linearly, 
while we sample parameters logarithmically if the variation exceeds this range. We provide 
an additional constraint on the parameters for the dark matter density profile, specifying 
that the total mass enclosed within 100 kpc fall between 0.25 and 2.0 times the value for 
the default Milky Way Halo. This is consistent with the fact that we select galaxies in our 
sample within a relativity narrow mass range as fixed by our constraint on Vfot- 

Figure 2 shows the distribution in the total luminosity calculated at 1.49 GHz for our 
Monte Carlo sample, illustrating the luminosity due both to dark matter annihilation as well 
as the cosmic ray luminosity for a galaxy with a cosmic ray injection rate equal to that of the 
Milky Way. A time-independent change in the cosmic ray injection rate would simply linearly 
shift the cosmic ray induced radio luminosity of a given galaxy. We also show for comparison 
the distribution of luminosities for all 66 "Milky- Way-like" candidates taken from the CA 
catalog, using the best-fitted distance given by NED. In Figure 3, we show the distribution 
of dark matter-induced luminosities as determined by our Monte Carlo as evaluated along 
the Galactic plane at a distance of r = 5 kpc from the GC (red line), and above the Galactic 
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Figure 2. Distribution of 1.49 GHz luminosity ratio L/ Lmw.dm where L is the total integrated 
luminosity and Lmw,dm is the luminosity due to dark matter for the canonical Milky Way model. 
The benchmark luminosity for dark matter plus cosmic rays is indicated by the vertical line. The 
simulated sample contains 2000 runs randomly distributed in the parameter space of tab. 3 with mass 
restricted to he within 0.25 — 2 times the mass of the Milky Way. We plot the contributions due to dark 
matter only and to cosmic rays only against extrapolated luminosities of 66 Condon Atlas galaxies 
meeting morphological cuts described in Section 2. Only 11.2% of the simulated dark matter- induced 
haze luminosities are a factor 20 smaller than our benchmark value, and only 1.42% are a factor 100 
smaller. 



plane at a height z = 1 kpc (blue line). Also repeated for comparison is the integrated 
luminosity due to dark matter. 

At this point, we note several key results: 

1. The total synchrotron luminosity from dark matter is almost invariably found to be 
subdominant to that from astrophysical cosmic ray sources, typically by more than one 
order of magnitude. We note from Figure 1 that this conclusion depends primarily 
on the dark matter profile, as the cosmic ray injection intensity is not altered in our 
simulations. This also implies that setting constraints on dark matter annihilation in 
external galaxies by assuming that the radio luminosity due to astrophysical sources 
is negligible is an extremely conservative approach, as the radio luminosities of most 
galaxies is expected to be dominated by cosmic rays from astrophysical sources. 
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Figure 3. Distribution of 1.49 GHz luminosity ratio L/Lmw.dm where L is the luminosity evaluated 
at a radius r — 5 kpc along the Galactic plane (red), height z = 1 kpc above the Galactic Center 
(blue) and Lmw,dm is the luminosity due to dark matter for the canonical Milky Way model. The 
sample is identical to that of Figure 2 and we have repeated the integrated dark matter component 
(green). 



4. 



Our modeled astrophysical synchrotron radiation has a more pronounced high lumi- 
nosity tail than is seen in the distribution of the CA galaxies, implying either that the 
range in the parameter space for cosmic ray diffusion is tilted towards very radio bright 
systems, or that the cosmic ray injection sources in our model galaxies are typically 
smaller than in the Milky Way. This is not true for our distribution of dark-matter 
induced emission, however, which is symmetrically distributed in log-space around the 
luminosity of the Milky Way. 

Only 11.2% of the simulated integrated dark matter- induced haze luminosities are a 
factor of 20 smaller than our benchmark value, and only 1.4% are a factor of 100 
smaller. This implies that the dark matter haze is a reasonably resilient feature among 
a large array of dark matter density profiles and diffusion scenarios. 

The simulated dark matter hazes can also be parameterized by the luminosity evaluated 
at 5 kpc along the Galactic plane and 1 kpc above the Galactic plane. In the r (z) 
directions, only 15.2% (12.3%) of simulations are suppressed by a factor 20 and 3.8% 
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(1.7%) are suppressed by a factor 100. In the next section we find that the strongest 
constraints are typically set by comparing the flux along the major- axis. For high 
inclination sources this is coincident with the radial direction considered here. 

5 Constraints from Radio Observations 
5.1 Methodology 

In addition to the above estimates of the total synchrotron power for a wide sample of 
values for the nuisance parameters, the Galprop code predicts in detail the morphology of 
the synchrotron emission stemming from dark matter annihilation. This allows us to set 
constraints on each CA candidate galaxy in two additional ways besides considering only the 
total integrated luminosity: 

• First, we can set constraints by comparing the peak flux contour in each galaxy against 
the predicted dark matter flux at the Galactic Center. 

• Second, we can compare the flux above or below the Galactic plane, where the dark 
matter produced haze is expected to be most significant compared to the synchrotron 
flux from cosmic rays. In order to do this, we calculate the fiux at the outermost 
observed radio contour for each galaxy and compare it to the modeled dark matter 
flux. We note that Figure 3 shows that the scale of the luminosity distributions from 
our dark matter models is expected to be similar in all cases. 

For the simulations presented in this section, we set all relevant parameters to the central 
values quoted in Table 3, except for the radial size size of the diffusion region, which we 
set to the size of the major axis as given in Table 2. We note that the diff^usion height 
employed in this analysis is extremely large, meaning that our choice of diffusion region is 
very conservative compared to the likely physical size of each galaxy. We then use Galprop to 
calculate the synchrotron emission throughout our model galaxy as a function of the radius 
and height, and integrate the emission over the line-of-sight based on the observed inclination 
of each candidate galaxy, correcting for the instrument beam-width for each galaxy. 

In order to compare the resulting flux contour map to observations of each CA candidate 
galaxy, we take into account the distance to each candidate galaxy and the related observa- 
tional uncertainties. Here, we only consider the central value for the distance measurement 
when the flux error due to the distance uncertainty is less than 25% of the standard devia- 
tion due to the systematic uncertainties discussed in Section 3. This condition does not hold, 
however, for the galaxies NGC 2683, NGC 4698 and NGC 4448. For these three galaxies, 
we provide flux contour-maps for the average, as well as the itlcr distance measurements. 
We note that in the case of NGC 4448, only three distance measurements exist in the litera- 
ture: we thus show results for the minimum, average, and maximum distance measurements 
reported in NED. For each model galaxy, we match our minimum contour to the minimum 
contour which is observationally reported, and then plot contours in flux increments of 2"/^ 
for integers n, as given by the CA catalog. For most galaxies, the lowest contour corresponds 
to n = —3, since the minimum observable flux/beam is relatively distance independent for 
resolved galaxies. However, the fluxes for Andromeda are reported with fluxes starting at 
n = -2. 

In order to calculate the peak flux contour from our Galprop models, we must employ 
a minimum angular region in which to calculate the peak flux, and this region must exceed 
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the size of the contour region resolved by the radio observations, or else we will set overly 
stringent constraints based on the smaller angular region considered in the Galprop models. 
We choose the peak flux to stem from the highest flux region with a physical size of at least 
25 kpc^, which exceeds the size of the radio beam at our maximum distance for all galaxies. 
We then round the reported contour down to the next contour level, using the 2"/^ contour 
spacing with the minimum contour flux set by the observation. We note that this assumption 
is conservative, as the peak radio flux within an observed beam is always slightly higher than 
the average beam over the matching Galprop contour. 

In calculations of the minimum flux contour along the minor axis of each galaxy, the 
anisotropy of the observed emission contours presents a challenge for any determination of 
the size of the candidate galaxies' minor axis. In order to conservatively estimate the size of 
the minor axis, we first calculate the direction orthogonal to the observed Galactic plane of 
each galaxy, and then calculate the smallest minor axis which includes all contours within 
10° of the minor axis. We note that in the case of NGC 2683 and NGC 7814, secondary 
sources are observed in the -z direction, and thus only the +z measurement is used to set 
contour constraints. This method effectively chooses the largest possible minor axis which is 
consistent with the flux data. 

In order to verify that our dark matter models are being correctly overlaid onto the 
radio halo, we also overlay the predicted dark matter contours onto optical images obtained 
by the STScI Digitized Sky Survey in Figure 4. In this case, we consider contours only for the 
central values of the distance to each candidate galaxy. We observe a strong morphological 
similarity, which serves to confirm our choices for both the size and orientation of the diffusion 
region. 

5.2 Constraints from the Integrated Flux 

In Table 4 we show the distance, integrated flux, and the estimate for the dark matter- 
induced radio haze flux as calculated using Galprop for our central parameter values. Also 
given is the ratio between the observed and simulated flux between 1.49-15 Ghz (Rmt) for 
each candidate galaxy. Ratios smaller than one, which we highlight in red, indicate that 
the observed flux is smaller than what we predict from dark matter annihilation alone. For 
galaxies with non-negligible distance uncertainties, the uncertainty in the distance, Galprop 
flux, and Rint are also given. 

We note two important results: First, Rint falls below unity for at least 15 of the 18 
galaxy observations, implying that the expected dark matter flux overproduces the entire 
observed radio flux from each galaxy, which presumably includes both dark matter and 
astrophysical sources of radio emission. While a constraint on dark matter annihilation from 
this simple process is difficult to quantitatively determine, we note that from Figure 2 it 
is clear that astrophysical cosmic rays usually dominate over the dark matter contribution, 
implying that for these candidate galaxies, drawn from a larger distribution of 66 Milky- Way 
like galaxies, the cosmic ray contribution is extremely low. Secondly, we note that for the 
only two galaxies with detected (non upper-limit) observations at multiple wavelengths (M31 
and NGC 2683), Rint decreases at higher frequencies in each case, demonstrating that the 
astrophysical emission may be significantly softer than the observed WMAP haze. 

Finally, we find that for five galaxies in our sample, Rint at 1.49 GHz is smaller than 
0.05. In our simulations of the dark matter haze, we find that obtaining a dark matter haze 
a factor 20 smaller than the benchmark values (i.e. i?int smaller than 0.05) has a probability 
of 11.2%. The likelihood of having at least 5 galaxies out of a sample of 66 with an emission 
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Figure 4. Overlays of dark matter only contours from Galprop onto optical images obtained from the 
STScI Digitized Sky Survey. Units are mJy beam~^. Contours are logarithmically spaced at levels 
2n/2 Jqj. jjjteger values n with outermost contours n — —2 for M31 and n — —3 for all others. With 
the exception of NGC 4448, which has very large distance uncertainties, the optical dimensions are 
similar or slightly exceed those of the lowest detectable radio contours. The diffusion zone has been 
matched to the estimated physical axis size. 

smaller than 0.05 corresponds to 87.5%, so we cannot derive any conclusion on the dark 
matter origin of the haze out of the integrated flux, barring stronger assumptions about the 
cosmic ray component of the radio emission. 
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Table 4. Integrated flux limits for selected galaxies at available radio frequencies from 1.49 — 15 
GHz. First and second values in parentheses represent limits corresponding to the minimum and 
maximum distance variations respectively. For measurements with associated flux uncertainties, i?int 
was computed using the 95% confidence level upper limit. Observations providing constraints are 
highlighted in red. 
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.125 


.0085 


.0043 


NGC1350 
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.006 
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.0054 


NGC2683 
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NGC4394 
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.0265 (.0162, .107) 


.0034 (.0033, .0034) 
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17.2 
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Table 5. Contour limits for peak, major (r), and minor (z) fluxes for selected galaxies at 1.49 
GHz. The first and second values in parentheses represent limits corresponding to the minimum and 
maximum distance variations respectively. Observations providing constraints are highlighted in red. 

5.3 Constraints from Peak and Axis Fluxes 

While we were not able to obtain any constraints from the total integrated radio luminosity, 
our dark matter radio emission simulations also allow us to investigate the signal for the 
galaxies in our sample by examining the peak emission (in the center of the galaxy) and 
off-peak emission (along and above the galactic plane). In Table 5 we show the flux ratios at 
1.49 GHz between each observation and our models at the location of the peak (Rpeak)) the 
minimum observed contour within 10° of the major axis (Rr), and the minimum observed 
contour on within 10° the minor axis (Rz). We note four interesting results: 

1. Evaluating the emission at specific locations, instead of the integrated emission, yields 
stronger constraints in every candidate galaxy. 
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Figure 5. 10 Kpc square overlay for M31 of Condon Atlas 1.485 GHz radio contours displayed over 
dark matter only predictions from Galprop. Units are mJy beam~^ and the image has been corrected 
for distance, inclination, position angle, and beam solid angle. The colormap is linear while contours 
are logarithmically spaced at levels 2"/^ for integer values n with outermost contours at the same 
level n = —2. Contour limits are computed using the conservative case (outermost observed contour) 
lying within 10° major and minor axis (r, z). Although the integrated and peak flux measurements 
do set limits for M31, the contour limits are the most stringent, constraining a dark matter haze at 
the < 1% level, relative to that of the Milky Way. 

2. The flux from dark matter annihilation in the majority of our candidate galaxies exceeds 
the observed radio flux (i.e. i?int < 1) at 1.49 GHz. 

3. We find that we get somewhat stronger constraints from R,. than from R^, somewhat 
contrary to expectations. In the case of the Milky Way the eccentricity in the ellipse 
stemming from dark matter annihilation is generally smaller than that for modeled 
astrophysical emission, so one might expect better constraints from off-plane observa- 
tions. 

In Figures 5 and 6 we show the observed 1.485 GHz radio overlays for M31 and for the 
remainder of our galaxy sample respectively, noting that the flux contours for each follow 
the form 2"/^ mJy beam~^ with an outermost contour of n = -2 for M31, and n = -3 for 
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Table 6. Fractions of random scans with luminosity suppressed by more than a factor of 20 and 
100 for a randomized sample of parameter space containing 2000 galaxies with masses restricted 
to .25 < M/Mmw < 2. Subscript "int" indicates the total integrated dark matter luminosity and 
subscript r (z) indicates the luminosity evaluated at a given distance along (above) the Galactic plane. 

the remainder of our candidate galaxies. In Figure 6 we show three distance variations 
corresponding to the minimum and maximum NED distance measurements for NGC 4448 
(2nd row) and the distance estimates corresponding to the la distance error for NGC 2683 
(top) and NGC 4698 (3rd row). The other three galaxies show a variation of less than 25% 
around the mean NED distance, and the variation is not considered. 

We note that in many cases the observed flux-contours are highly anisotropic, which 
contrasts with the approximately symmetric contours expected from dark matter. While 
anisotropics could be introduced by the cosmic ray propagation parameters, the one-dimensional 
simulations from Figure 1 show that these anisotropies in the magnetic field parameters and 
diffusion constants have relatively minor effects on the radio luminosity, compared to the 
case of astrophysical injection. Instead, the odd shapes of several galaxies' contours, partic- 
ularly those of NGC 4698 and NGC 7814, are highly indicative of anisotropic astrophysical 
cosmic ray injection. Most importantly, we do not see elliptical emission profiles for any 
galaxy, contrasting with our model for the Milky Way, which has a significant injection of 
cosmic rays in star formation regions across the Galactic plane. This means that in these 
radio-underluminous candidate galaxies, the elliptical extension of the expected dark matter 
profile in the radial direction produces significantly better limits at than in the direction 
of the minor axis. 

The results of Table 5 indicate that for three galaxies the off-peak emission along the 
galaxies' major axis (r direction) is suppressed by more than a factor 100 with respect to our 
benchmark dark matter haze model. Our simulations, reported as the red line in Figure 3 
and in Table 6, indicate that only 2.3% of the simulated galaxies are suppressed by at least 
.006 compared with the benchmark dark matter model. We calculate that having at least 3 
galaxies out of the 66 candidates with such a suppressed radio emission has a likelihood of 
less than 20%. Although it is difficult to turn this estimate into a confidence level exclusion 
limit to the dark matter haze hypothesis, the large variations we allow in multiple physical 
parameters make it clear that there is some tension in the dark matter interpretation of the 
galactic haze with observations of external spiral galaxies. 

5.4 Future directions 

With upcoming data from the Planck satellite, we expect to obtain a much more robust 
analysis of the Milky Way's haze. Planck will produce all-sky observations over nine frequency 
bands from 30 to 857 GHz, providing a much needed spectrum with superior foreground 
identification and separation, as well as detailed polarization data. Having established an 
effective lower bound on the integrated dark matter luminosity, the availability of an all- 
sky spectrum up to higher frequencies will allow us to transcend the issues associated with 
complicated cosmic ray sources. With an angular resolution of 30 arcminutes at 30 GHz and 
24 arcminutes at 44 GHz, Planck will have the ability to resolve galaxies up to roughly 5 Mpc. 
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Figure 6. 20 kpc-square overlays of Condon Atlas 1.485 GHz radio contours onto dark matter only 
predictions from GALPROP. Units are mJy beam~^ and images have been individually corrected 
for distance, inclination, position angle, and beam solid angle. For NGC 4448, distance variations 
correspond to minimum and maximum NED measurements centered on the value from the Condon 
Atlas. For NGC 2683 and 4698, distances are allowed to vary by 1 standard deviation from the mean 
based on NED measurements. The colormap is linear while contours are logarithmically spaced at 
levels 2"/^ for integer values n with outermost contours at the same level n = —3. A dark matter 
haze should produce a fairly symmetric synchrotron signal. Here the irregular shapes of the observed 
contours often yield significant constraints along r, z, or both. Additionally, the peak predicted 
contours are often substantially higher than what is observed. 
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This encompasses the local group which contains several well studied Milky Way analogs. 
Galaxies at larger distances will still register as point sources. If we exclude sources with 
known nearby radio confusion and cross-reference against a larger galactic catalog to make 
morphological cuts, it should be possible to obtain integrated flux limits with substantially 
better statistics than those presented here. 

In order to test the impact of Planck observations, Figure 7 shows the predicted total 
luminosity components at 30 and 44 GHz, compared to the benchmark Milky Way value 
for galaxies in our random parameter scan. The distribution is similar in width to that 
of the 1.49 GHz distribution shown in Figure 2. However, the synchrotron spectrum from 
dark matter is expected to be significantly harder than that of cosmic rays. Thus, at higher 
frequencies, we expect a better dark matter signal-to-cosmic ray noise ratio. Additionally, 
at 30 GHz (44 GHz), 6.49% (1.4%) of models produce a haze with more than a factor 100 
suppression compared to the Milky Way model. With a larger sample, higher quality data, 
and higher-frequency observations, it should be possible to definitively rule out dark matter 
annihilation as the primary progenitor of the haze, even under very conservative assumptions. 

6 Discussion and Conclusions 

If the WMAP haze is indeed produced by synchrotron emission from electrons and positrons 
generated via dark matter annihilation, galaxies with properties similar to those of the Milky 
Way should also possess an observable radio haze. High-resolution radio data allows for mor- 
phological studies in addition to peak and integrated flux comparisons against simulations. 
In this paper, we have examined seven galaxies of size and morphological type similar to 
the Milky Way, and compared the available data at 1.49-15 GHz against the predicted dark 
matter haze, as calculated using the code Galprop. Our results are conservative in the sense 
that we did not subtract the cosmic ray-induced component of the radio emission, which is 
typically dominant, when setting limits on the dark matter synchrotron emission from the 
galaxies under consideration. 

Although there are significant systematic uncertainties involved in this procedure, our 
analysis found that most of the galaxies under consideration are less luminous at radio fre- 
quencies than would be expected from their dark matter-induced haze, unless somewhat 
extreme parameters are adopted for their dark matter halo profiles, cosmic ray propagation 
models, or magnetic field models. In particular, to reduce the dark matter luminosity below 
the observed upper limits, one must simultaneously reduce the dark matter density, increase 
the diffusion constant, and/or increase the interstellar radiation energy density relative to 
those features of the Milky Way. The size of the diffusion region and the strength of the 
magnetic field are less important to this problem. While most of these parameters are very 
poorly constrained, we have performed a Monte Carlo scan over what we consider to be a 
reasonable range of the parameter space, and find that the dark matter-induced haze is less 
than 1% as luminous as the Milky Way's haze in only about ~ 5% of the simulated galaxies. 

To obtain limits on the dark matter origin of the haze, we compare observed and pre- 
dicted flux densities along the semi-major and semi-minor axes, peak contour levels, and 
integrated fiux. For six of the seven galaxies examined, a dark matter interpretation is con- 
strained below the 5% level. We obtain stronger constraints from observations at a position 
of r=5 kpc along the galactic plane, where we find three galaxies which are underluminous 
by a factor of 100 compared to the Milky Way. Qualitatively, it is also worth noting that 
the radio emission from the majority of these galaxies does not closely resemble the mor- 
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Figure 7. Distribution of 30 and 44 GHz luminosity ratios L/Lmw,dm where L is the total integrated 
luminosity and Lmw,dm is the luminosity due to dark matter for the canonical Milky Way model at 
the respective frequency. Sample contains 2000 runs randomly distributed in the parameter space of 
Table 3 with the mass restricted to lie within 0.25 to 2 times that of the Milky Way. We plot the 
contributions due to dark matter and cosmic rays individually. At both 30 and 44 GHz, 6.5% and 
1.4% of models show suppression greater than a factor of 100 compared to the benchmark Milky Way 
model respectively. This indicates that higher frequency observations from Planck will produce more 
robust limits on dark matter origin of the haze. 



phology expected from dark matter annihilation, and instead appears to trace spikes in the 
cosmic-ray injection rate. Remarkably, one of the galaxies without evidence of a dark mat- 
ter haze is our local-group neighbor, M31 (the Andromeda Galaxy). While it is difficult to 
turn this into a firm exclusion limit on the dark matter origin of the Milky Way haze, one 
would have expected Andromeda to be significantly more radio luminous in most dark matter 
interpretations of the WMAP haze. 

Upcoming data from the Planck space observatory will offer a much broader set of ra- 
dio frequency observations (nine bands from 30 to 857 GHz). Notably, recent work from the 
Planck collaboration has already confirmed presence of the WMAP haze in the Milky Way, 
and found that it has a microwave edge spatially coincident with the gamma-ray bubbles 
detected by Fermi. Such a broad frequency range, and additional polarization information, 
will allow for much greater separation of astrophysical foregrounds, as well as reduced con- 
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tamination from cosmic rays. 
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